Collective Modes in the Loop Current Ordered Phase of Cuprates 
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Recently two branches of weakly dispersive collective modes have been discovered in under-doped 
cuprates by inelastic neutron scattering. Polarization analysis reveals that the modes are magnetic 
excitations. They are only visible for temperatures below the transition temperature to a broken 
symmetry phase which was discovered earlier and their intensity increases as temperature is further 
decreased. The broken symmetry phase itself has symmetries consistent with ordering of orbital 
current loops within a unit-cell without breaking translational symmetry. In order to calculate the 
collective modes of such a state we add quantum terms to the Ashkin- Teller (AT) model with which 
the classical loop current order has been described. We derive that the mean field ground state of 
the quantum model is a product over all unit-cells of linear combination of the four possible classical 
configurations of the loop current order in each unit-cell. The collective modes are calculated by 
using a generalized Holstein-Primakoff boson representation of orbital moment operators and lead to 
three branches of gapped weakly dispersive collective modes. The experimental results are consistent 
with the two lower energy branches; the third mode is at a higher energy than looked for by present 
neutron scattering experiments and might also be over-damped. Implications of the discovery of the 
collective modes are discussed. 

PACS numbers: 



I. INTRODUCTION 

Thermodynamic 1 ^ as well as other properties 3 in all 
Cuprate families with high superconducting transition 
temperatures change their temperature dependence be- 
low about a characteristic temperature T*(x), which de- 
pends on doping a;, from those at higher temperatures. 
The region below T*(x), see Fig. (JIJ is said to mark 
the pseudogap region. Photoemission 4 " and Angle re- 
solved photoemission experiments reveal^ an anisotropic 
decrease in single particle spectral weight near the chem- 
ical potential below T*(x). A popular belief^ has been 
that T* (x) marks a crossover to a region of reduced low 
energy one-particle spectral weight, as well as in multiple 
particle-hole spectral weights, due to one or other rea- 
son, preformed pairs^, resonating valence bonds^, stripe 
formation 8 or other states of charge modulation, proxim- 
ity to the Mott-AFM-insulating state 9 , etc.. Guided by 
the fact that a "strange -metal " region whose proper- 
ties can be explained by quantum-critical fluctuations^ 
abuts the pseudo-gap region in the phase diagram, a 
broken symmetry state was sought for the region below 
T*(x) with T*(x) —¥ at x — > x c , the quantum-critical 
point. An unusual class of phase d 11 ' 12 was found to be 
stable in mean-field calculations of the three-band model 
for cuprate a 14 ' 15 . In such phases time-reversal symmetry 
is broken via spontaneous generation of orbital currents 
in each cu-o unit-cell without altering the translational 
symmetry. The phase transition belongs to the class with 
an order parameter singularity but no divergence in the 
specific heatiSii 2 ^. 

An order parameter consistent with the symmetry 
predicted has been discovered in four different families 
of cuprates through polarized neutron scattering^ or 
dichroic ARPES 1 -^ and may be regarded as a universal 
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Figure 1: The universal phase diagram of hole-doped cuprates 
based on properties which show characteristic changes across 
the boundaries shown in all the cuprates. The boundary of 
the pseudogap phase has not yet been determined in experi- 
ments, i.e how T*(x) line continues for lower x has not been 
determined. 



property of the cuprates. The magnitude of the order 
parameter is large, estimated in ff<?1201, with T c w Q1K 
and T* s» 350-rY, to be about OA/ib per unit-cell. 

More recently inelastic neutron scattering has 
discovered 1 ^ two branches of weakly dispersive collective 
modes for three different x in the family (Hgl20l) for T 
below T*(x). (Similar modes have been also found 1 ^ in 
Y BclzCuzOq.q.) Their onset temperature and increase of 
spectral weight as a function of temperature follows the 
temperature of the intensity of the elastic order. Such 
modes were also anticipated 1 ^ and further substantiate 
the broken symmetry. In this paper, we describe details 
of the calculation of such modes as well as of another 
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Figure 2: The four Possible "classical" domains of the loop ordered state are shown. In the classical ordered phase, one of these 
configurations is found in every unit-cell. 



higher energy branch of excitations which has not yet 
been discovered. A brief report of this work has already 
been published^. 

The observed broken symmetry is consistent with 
spontaneous moments due to a pair of orbital current 
loops within each unit-cell preserving overall transla- 
tional symmetry. It breaks both time-reversal and inver- 
sion symmetry, preserving their product. The "classical" 
order parameter^ may be characterized by the anapole 
vector— L 
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where the moment distribution M(r) is formed due to 
the currents on the four O-Cu-0 triangles per unit-cell 
as shown in Fig. @. This figure also shows the four 
possible "classical" domains of the loop current ordered 
state. In the classical ground state, ordering occurs in 
one of the domains shown. 

Quantum-mechanics allows local fluctuations among 
the four configurations in Fig @. This leads, as shown 
in this paper to a ground state in which each unit-cell has 
a finite admixture of all the four configurations. It also 
leads to three branches of collective modes of the order 
parameter at finite energies at all momenta q for T < T* . 
The finite energy follows from the fact that the ground 
state has symmetry consistent with that of a generalized 
(transverse-field) Ising model. In this paper these modes 
will be derived. One can argue that there should be three 
because each of the four configurations can make transi- 
tions to the other three as pictorially shown in Fig. ([3]). 

This paper is organized as follows: In the next sec- 
tion, we introduce the classical AT model for the loop 
current order and generalize it to the quantum model in 
the SU(4) representation rather than the SU{2) x SU(2) 
of the classical AT model. The quantum terms are cho- 
sen from considerations of the internal and lattice sym- 
metries of the classical model. In the following section, 
the ground state of the quantum model is evaluated in 
mean-field and the dispersion is calculated using the gen- 
eralization of the Holstein-Primakoff transformation. We 
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Figure 3: The schematic figure shows that there are only 3 
collective modes. 



compare with the results from experiments. We conclude 
by discussing the significance of the experimental discov- 
ery of the collective modes and the further possible effects 
which arise from the calculations here. In four Appen- 
dices, we discuss the necessity for casting the problem in 
the SU (4) representation, some technical details, and the 
theory for inelastic neutron scattering from the collective 
modes. 



II. MODEL FOR QUANTUM-STATISTICAL 
MECHANICS OF LOOP-CURRENTS 



The order parameter L and an effective Hamiltonian 
for this collective variable has been derive d 11 ' 12 ^ 3 start- 



ing from a model of interacting fermions^. The full 
Hamiltonian may be written as a sum of three parts: 



H = H r 



Hp 



(2) 



where H co u is the Hamiltonian for the collective co- 
ordinates which order at T* to give the long range 
order L, Hp is the remnant fermion Hamiltonian and 
Hp-coll i s the residual interaction between the fermions 
and the collective coordinates. In this paper we will be 
concerned almost exclusively with H co u, although in a 
brief discussion of damping near the end, Hp- co u is im- 
plicated. 

An approximate representation of L of Eq. (JlJ is given 
by M(r) « M M c5 2 (r - r M ), where r M , = 1,..,4 are 
the location of the four "sites" in any cell at the cen- 
troid of the moment distribution. These sites are labeled 
S\,..,Si in Fig. (01 . The orbital moments are ei- 
ther up or down or zero. The four classical domains of 
Fig. may be represented by the four values of the angle 
9 = 7r/4, 37r/4, 5tt/A, 7tt/A that L makes with the x— axis. 
The 4 classical loop current states are eigenvalues of op- 
erators L* — {L itXl Li^ y ) defined at the unit-cell i. We 
may define a basis, choosing |Lj| to be unity, 



(L i x + iLi v )\d)i = e 



(3) 



The classical statistical mechanics of the Loop-Current 
state may be derived from the Ashkin- Teller model, 
which is given in terms of a pair of Ising spin per unit-cell 
of,T?, whose eigenvalues, ±1 specify the x and y com- 
ponents of the direction of the vector L. The four loop 
current states can therefore also be denoted as | ± 1, ±1). 
The classical Ashkin- Teller model 2 ^ is given by2£, 

H AT = ~ [«^i <7 f cr i + ^rfrf + ^^ajr/r/] (4) 
(i,j) 

Quantum fluctuations among the four possible direc- 
tions of order together with dissipation lead 26 to a scale 
invariant spectrum which leads to the observed Marginal 
Fermi-liquid properties in the quantum-critical regime in 
the phase diagram, Fig. (JXJ) . In this paper, we will derive 
the effect of the quantum fluctuation in the ordered loop 
current phase. In this phase, the kinetic energy locally 
flips a loop current state in cell i to one of the other three 
states, just as the transverse field does between the two 
states of the transverse field Ising model. Due to this 
term, the ground state is a superposition, in each unit- 
cell, of all the 4 possible directions of Lj. The excitations 
consist of local flips between the configurations which 
spread out spatially through interactions between neigh- 
boring sites and acquire dispersion just like spin waves 
in the transverse-field Ising model. 



A. Symmetries 

As explained in Appendix B, in order to treat the four 
states on the same footing in the quantum-mechanical 




Figure 4: We represent here the notation introduced in the 
text for the location of the moments within a unit-cell. 



model, we must enlarge the representation from that of 2 
by 2 matrix space to that of 4 by 4 matrices, i.e. we must 
consider the problem in the ST/ (4) representation rather 
than in the SU{2) x SU(2) representation. We introduce 
the following direct products between Pauli matrices. 



S 4 = a 1 <g> I, T l = I®t\ K lJ =(7*0 r j 



(5) 



The 15 traceless matrices S, T and K are also the gener- 
ators of SU (4) algebra. Their commutators are easy to 
compute. 

[S\ S J ] = 2ie ijk S k , [T\ T J ] = 2ie ijk T k , (6) 
[S\TJ]=0, ( 7 ) 
[S\ K jk ] = 2ie tjl K lk , [T\ K^ k ] = 2ie ikl K^ 1 , (8) 
[K*S, K kl ] = 2ie ikm S m 5 jl + 2ie jln T n 6 ik (9) 

Under this choice of basis, the 4 states are labeled by 
the eigenvalues of a z and r 2 . We will use the following 
short-hand notation 



1,1) 



ID 



|l>r 



(1,0,0,0) T . (10) 



And similarly for the other 3 states | — 1, 1), |1, — 1), | — 
1,-1). In this basis, the classical AT model can be 
rewritten as 

H = -J2 l J iSfS* + J 2 TfT* + J A KfKf] (11) 

(i,3) 

We now list both the internal and the lattice symmetries 
of this classical model. 



1. Internal symmetries: 

The classical AT model has an obvious global contin- 
uous symmetry 1/(1)53 x J7(l)y3 x U(1)k 33 which is the 
rotation around S 3 , T 3 and K 33 generated by the rota- 
tion matrix U(0) = e^ 5 ' / 2 e ^ T3 /V^ 33 /a _ xhis is the 
only continuous symmetry of the model and implies that 
the phase differences between the four states are arbi- 
trary. Any quantum term will break this symmetry. 
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The classical model also possess several discrete sym- 
metries. We will argue that the quantum terms should 
preserve the discrete symmetries. 

Like the Ising model, the classical AT model has a 
global Z2 symmetry which maps tr 2 to — a z . Since all 
the interaction terms involve two a z , they stay invariant. 
Similarly, we also have anther Z2 symmetry which maps 

T Z tO — T Z . 

We are more interested in the symmetric case, J\ = 
J2. In this case, there is also a discrete symmetry Z4 
which is the symmetry under cyclic permutation among 
the 4 states. Since the coefficient of S 3 and T 3 terms 
are the same, the Hamiltonian is the invariant under this 
transformation. There is another Z2 symmetry which 
switch a and r. 



as the internal Z2 symmetry a r and is generated by 

/ 1 \ 
-/ I 1 I ,,,, 

°W4) = oiOO ^ 

\ 1 / 

Similarly, the reflection about y — —x switches |1, 1) and 
I — 1,-1) and keeps the other two states the same. It is 
generated by <r(37r/4) = S 1 T 1 <t(tt /4) which is a combi- 
nation of the 3 internal symmetries. One can see that all 
the internal discrete symmetries are coincident with the 
lattice symmetries. We expect that the quantum term 
preserve all the lattice symmetries, so it also preserve all 
the internal discrete symmetries. 



3. Quantum terms allowed. 



2. Lattice symmetry 

The 4 states represent the loop current pattern on the 
square lattice of copper and oxygens. The point group 
of a square lattice is D4. It consists of 4- fold rotations 
and reflection about the x-axes, the y-axes and the two 
diagonal axes. The loop current states can be thought 
as currents flow along both x and y axes and the values 
of a z and t z label the directions of these currents. From 
these, one can find out how the point group of lattice act 
on these 4 state. 

The 4-fold rotation makes the following transformation 
|1, 1) ->■ I - 1, 1) -+ I - 1, -1) ->• |1, -1) ->■ |1, 1) which is 
the same as the internal Z 4 symmetry. The transforma- 
tion matrix is given by 



C 4 = 



/o 1 
10 
1 

Vo 1 



(12) 



If we think of |1,1) as a plane vector (1,1), then the 
above transformation makes a n/2 anti-clockwise rota- 
tion, which realizes the cyclic permutation among the 4 
vectors. Then one can find out the following transforma- 
tion 



Now we want to identify possible terms which lead to 
the quantum fluctuations among the 4 classical states. 
Ir>2&, the AT model has been mapped to XY model in 
the quantum-critical fluctuation regime of the phase di- 
agram. In the loop ordered state, quantum term of AT 
model can be thought as a discrete version of the quan- 
tum term of XY model which makes clockwise and anti- 
clockwise rotations. As one can see, the unitary oper- 
ator to generate anti-clockwise rotations is just C4 we 
introduced. But since detC4 — — 1, this operator is not 
a proper rotation in the complex 4D space but also in- 
cludes a parity transformation e.g. |1, 1) — > — 11, 1) and 
keep other 3 states which cannot be realized as a con- 
tinuous transformation. Since the quantum evolution of 
states is obtained by integrating infinitesimal rotations, 
we expect the anti-clockwise rotation operator U to be a 
proper rotation, or dctU = 1. (This matter is further dis- 
cussed in Appendix A). Therefore, we introduce a phase 
factor to correct the sign of the determinant. 



U = 



/o 1 
1 1 W4 

1 1 ' 

Vo 1 



(15) 



Then the possible quantum terms should be all the pos- 
sible Hcrmitian combinations of this unitary operator: 



C 4 S 3 C] = -T 3 



CaT 3 C\ = S 3 



CaK 33 C\ 



-K 3 fl3) 



[/ + t/ f 



1 



V2 



[S 1 + T 1 + (K 26 - K i2 )) (16) 



Clearly, the classical AT model is invariant under this 
transformation. 

Now we turn to the reflections. The reflection about 
x-axes switches | ± 1,1) with | ± 1,-1), which is also 
equivalent to the Ising like internal Z2 symmetry. This 
reflection can be generated by operator cr(0) = T 1 Sim- 
ilarly, the reflection about y-axes switches |1,±1) with 
I — 1,±1) and is generated by operator a{ir/2) = S 1 . 
The reflection about y — x switches |1, — 1) and | — 1,1) 
and keep the other two states the same. This is the same 



i(U-rf) = -j=[S 1 -\ 
l -[U 2 -{U^) 2 ]=K^ 



- T 1 - (K 23 - K 32 )] (17) 
(18) 

Since U 4 = -1, we have U 2 + (W) 2 = 0. Therefore the 
above are the only possible independent Hermitian com- 
binations to quadratic order in U. In the most general 
case, all these terms should be present in the quantum 
AT model. For our problem, we need to constrain this 
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further to be consistent with the lattice symmetry. The 
4-fold rotation operator C4 commutes with U, thus there 
is no constraint due to 4-fold rotation. As discussed, the 
reflection operators <r(0) and a(ir/2) are equivalent to T 1 
and S 1 , respectively. And, for the reflection about x and 
y axes and diagonal directions, we have 

S 1 (U + U r )S 1 =i(U- 17+) (19) 
T^U + U^T 1 = i(U- U^). (20) 

For reflection operator cr(±7r/4), we note that 

a{n/A){U + U^)d{n/A) = i(U-U*) (21) 
ct(3tt/4)(?7 + U*)a(3ir/4) = i(U - U f ). (22) 

K 11 is invariant under all these 4 reflections. Once 
can see that S 1 + T 1 is even under the reflections and 
K 23 — K 32 is odd. Therefore, the lattice symmetry re- 
quire the coefficient in front of U + TP and i(U — TP) 
to be the same in order to cancel out the odd part. In 
conclusion, the proper quantum term that respect all the 
lattice symmetry should be S 1 + T 1 and K 11 . 

Even if we do not start with the rotation operator U, 
we will get the same results for the quantum term just by 
imposing lattice symmetry. In the chosen basis, the only 
diagonal matrices are S 3 , T 3 and K 33 and they form a 
Cartan sub-algebra^ of the SU (4) Lie algebra, which is 
the maximal commuting sub-algebra one can find. All 
the other 12 generators contains non-zero off-diagonal el- 
ements and act as raising and lowering operators, which 
mix among the 4 states. The most general quantum term 
should be the superposition of all these 12 generators. 
Lattice symmetry requires that the quantum term com- 
mute with S 1 , T 1 , K 11 and <t(tt/A). It is easy to see 
that S 1 , T 1 , K 11 commute with each other thus they 
also form a Cartan sub-algebra. Therefore, the only gen- 
erators which commute with S , T , K 11 are themselves. 
ct(7t/4) switches S 1 with T 1 , thus require the coefficients 
of S 1 and T 1 to be the same. Thus the quantum term 
S 1 + T 1 and K 11 is unique up to some unitary transfor- 
mations. 

It may seem strange that only the generators with 1 
component are involved in the quantum term. Actually, 
this is only due to the specific choice of representation 
we have made. We can redefine a(0) and a(w/2) as 
ct(0)|1,1) = e^|l,-l) and ct(tt/2)|1, 1) = e+*\ - 1,1) 

a(0) =I®( e \ e ^ J = T 1 cos 0! - T 2 sin ^(23) 

a(7r/2) = ^ e \ J ® I = S 1 cos<f> 2 - S 2 sin0 2 .(24) 

Since the D4 group can be generated by 2 elements, once 
the phase factor in <r(0) and a{n/2) are fixed, phase of all 
other D4 operators are also determined. In general, there 
are three arbitrary phase difference one can choose, but 
in the symmetric case, there are only two left. The other 
symmetry operators are obtained from C 2 — a(w/2)a(0), 
ct(tt/4) = C 4 ct(0) and a(37r/2) = C*|ct(7t/4). As before, 



cr(0), a(w/2) and cr(7r/2)cr(0) form maximal commuting 
sub-algebra. Therefore the quantum terms are given by 
ct(0) + a(n/2) and a(n/2)a(0) with both x and y com- 
ponent present. The new cr(0), a(n/2) are related to the 
old ones by unitary transformations 

e t(4>iT 3 +4, 2 s 3 )/2 s i e -i(<t >1 T 3 +<t >2 S 3 )/2 = o-(7r/2) 

e i(0 1 T 3 +0 2 S 3 )/2 T l e - 4 (0 1 T 3 +0 2 S 3 )/2 = ~/qS ^5) 

Thus the results for energy of the collective modes (but 
not their eigenvectors or the ground state wave- function) 
should remain the same as in the original gauge. 

III. DISPERSION OF COLLECTIVE MODES 

In the SU(4) formalism, the classical Ashkin Teller 
model can be written as 

H AT = - [JiSfS 3 + J 2 T 3 T 3 + J 4 K 33 K 33 } (26) 

According to the discussion of section [III the quantum 
fluctuations can be described by the Ashkin Teller model 
with a symmetric quantum terms or transverse field as 
S 1 + T 1 and K n . The Hamiltonian for the quantum 
Ashkin- Teller model is therefore given by 

i 

- ^[JiSfS? + J 2 T 3 T 3 + J A K 33 K 33 } (27) 

As discussed earlier, the basis of 4 states at each lattice 
site correspond to the 4 types of loop current states. In 
the AT model, these 4 states are labeled by the eigen- 
values of a z and t z which can be denoted as | ± 1, ±1). 
The kinetic term S 1 flips the sign of first index and T 1 
flips the second. Then K 11 flips both indices together. If 
we think of the two indices as x and y components of a 
plane vectors, the 4 states can also be represented by 4 
plane vectors. Then S* 1 and T 1 rotates the vector by ir/2 
and K 11 rotates the vectors by ir. In principle the pa- 
rameters in this model should be determined from micro- 
scopic models. We will treat them as phenomenological 
parameters and determine them by fitting the observed 
dispersion of the collective modes, and judge if they have 
the scale of values expected, which is of the order of the 
pseudo-gap energy. 

A. The Ground State 

In order to calculate the collective mode dispersion, 
we first determine the ground state in the ordered phase 
at the mean field level. For the loop current ordered 
phase the classical AT model has ferromagnetic interac- 
tions J1.2 > and J4 < |Ji,2|, the latter to ensure that 
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there is no divergence of the specific heat at the tran- 
sition. In this phase, we have (S 3 ) 7^ 0, (T 3 ) ^ and 
(K 33 ) 7^ 0. The most general uniform quantum ground 
state is a product over all sites of the linear combination 
of the 4 basis states. Since the overall complex factor 
is redundant, one can parameterize this state by 6 real 
parameters. Thus we assume the ground state to be 

|G)o = p| cos ^- ( cos ^|l, 1), + sin ^ |1, -1>,) 

i 

+ sin|e*(cos^|-l,l). i 

+ sin | e *|_ 1) _ 1) .) (28) 

It is straightforward to compute the energy per site E g 
but the result is very complicated and we have to use 
numerics to minimize E g to find out parameters 9i and fa . 
Since the Hamiltonian does not involve the 2 component 
matrices and only 2 component matrices are imaginary, 
we expect the 3 phase parameters fa — for i = 1,2,3. 
Thus all the 4 coefficients are real numbers. This can 
be checked by numerical minimization of the energy. In 
comparing with experiments, we are more interested in 
a special case that J\ — J2 = J and the Hamiltonian is 
invariant under the interchanging between tr and t. In 
this case, we still have to numerically minimize E g , but 
we find that the result can be well approximated by a 
direct product state 

\g) = n( c °4 |i>i+si 4 |_i> 

j \ Jo 
® ^COB^|l)i + Bill||-l)(J (29) 

which is parameterized by one angle 0. Then the angle 
is determined by equation, 

t + t' sin 9 + 4 J sin 9 + 4 J 4 sin 9 cos 2 9 = 0. (30) 

We will use this direct product ground state |G) for lat- 
ter discussions. The energy of this direct product ground 
state is only 3% higher than the most general ground 
state, so it is a good approximation. We emphasize that 
this wave-function is in a particular gauge; the change 
of gauge, Eqs. (|2"5j) . changes the ground state wave- 
function. 



B. Collective Modes 

To compute the spin wave like collective modes in 
this quantum Ashkin- Teller model, we generalize the 
Holstein-Primakoff transformation for SU(2) to SU(4). 
Introducing three boson operators bi, Ci and di, we 
have the following boson operator representations for the 
Ashkin- Teller model (For ease of notation, we omit the 



lattice site label i ): 

S 1 = Uc + c)U + bU+dtb, 

S 3 = l- 2c f c - 2dU, 

T 1 = Ub + tfU + Jd + d?c, 

T 3 = 1 - 26 f 6 - 2dU, 

A 11 = Ud + cftU + tfc + Jb, 

A 33 = 1 - 26+6 - 2c t c, 

K 13 = Uc+JU -tfd-Sb, 

K 31 = Ub + tfU-cU-dtc, 

with U = (1 - b f b - c*c - d f d) 1/2 . In this repre- 
sentation, we take the classical ground state such that 
(S 3 ) = {Tf) = (K 33 ) = 1. Thus the classical ground 

state is just |1, 1) and b\, cj, d\ are the creating opera- 
tors of states |1, —1), I — 1, 1) and | — 1,-1) at site i from 
the classical ground state respectively. 

Due to the quantum terms, the mean field ground state 
\G) is a superposition of the 4 states at each site. In order 
to make use of the Holstein-Primakoff transformation, we 
can make a basis rotation to transform the ground state 
\G) to rj s |1, l)i in the new basis. 

In the following calculations, we still keep J\ and 
J2 as two different parameters. In the last step, we 
will take J\ = J2 = J. For the direct product 
ground state \G) = FJi ( cos IH 1 )* + sin |H ~ l )i) a ® 
(cos 4f + sm^-| — , the rotation we need is just 
two rotations around y axes by angle 9 1^2 in the a and r 
space. 

S 3 -> S 3 cos 0i - S 1 sin 0j T 3 -> T 3 cos 9 2 ~ T 1 sin 9 2 
S 1 -> S 3 sin 6»i + S 1 cos 9 X T 1 -> T 3 sin 9 2 + T 1 cos 9 2 

Then the Hamiltonian is 

H = + Sla) + t(T 3 s 2 + T/ca) 

i 

+t'(K 33 Sl S2 + K 31 Sl C 2 + K} 3 Cl S2 + K\ Y Cl C2)\ 

- J2 [ J i(S?S 3 cl - S 3 Sjc lSl - S\S)dS\ + S\S)s\) 

i,j 

+J 2 {T 3 T 3 cl - T 3 T]c 2 S2 ~ T}T 3 c 2 S2 + T}T]s\) 

-J^{K 33 c lC 2 - K 31 c lS2 - K 13 Sl c 2 + K 11 s 1 s 2 ) i 

id 

x (A' 33 cic 2 - A 31 cis 2 - A 13 sic 2 + # u sis 2 )j. (31) 

with d,2 = co s 0i,2 an d si,2 = sin.0i 2< The 0i ; 2 can be 
solved from the equations. 

t + t's 2 + 4JiSi + 4J4S1C2 = 
t + t's 1 + 4J 2 s 2 + AJ A s 2 c\ = 

1. Simple Case, J4 = t' = 

We can start from the simple case with J4 = and 
t' = 0. In this case, the Hamiltonian is two decoupled 
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Ising models. We can plug in the boson transformation, 
expand the square root up to quadratic order. Since we 
only care about the collective mode dispersion, the con- 
stant term can be ignored. The linear terms in boson 



operators will cancel out due to the minimization condi- 
tion of the ground state energy. Thus we are only left 
with quadratic terms, 



H =*^[(-2c t c-2d t d)iSi + [bU + Sb) lCl ] - J x ^ [(-2^0- 2dU)ic\ + (-2c f c- 2dU)jc\ 

-(tfd + <i t 6) i sici - (tfd + d^fyjSxC! + (c + c f )i(c + c^s 2 + terms (b O c, 1 -H- 2). (32) 



This equation can be simplified by using identity — tsi + 4Jicf — 4Ji. Note that (b^d + d^b) and (c^d + d^c) terms 
cancels out, thus operator d decoupled from b and c . In momentum space we have 



H = X]( 4J i ~ 2J i s i/k)(4 c k + c-kdk) - 2J lS 2/ k ( CkC _ k + c l c l k ) 

k 

+ (4J 2 - 2J 2S 2/ k )(6t 6k + fe_ k 6t_ k ) _ 2J 2 s 2 / k (6 k 6_ k + 6 k & f _ k ) + 8(Ji + J 2 )4^k. 



(33) 



One can see the that the three sets of 3 bosons are de- 
coupled. The b and c part are two copy of Ising models 
and the d part is free boson. Making use of standard 
Bogoliubov transformation we can diagonalize the above 
Hamiltonian and find the following dispersions 

^=2[(4J 1 ) 2 -t 2 / k ] 1/2 , 
^2 = 2[(4J 2 ) 2 -t 2 / k ] 1 /2, 

oj 3 = 8(Ji + J 2 ). 

In addition to the two Ising model modes, we also find a 
third mode which is dispersion-less. The two Ising modes 
are generated by bosons operators b and c which only 
rotate a or t space alone. Since the ground state is a 
direct product, the excited state generated by modes b 
and c are also direct products. 

exp(-#6)|l,l) « 

= |1> x (|1> - 1>). (34) 

On the other hand, the third mode which is generated by 
boson operator d cannot be written as a direct product. 
For example, 

exp(-tyd)|G ) w |1, 1) - i4\ - 1, -1), (35) 

which cannot be written in the form: \ip a ) <8> IVv)- 

The constant energy feature of the third mode can be 
explained as follows. The Hamiltonian with J4 = in 
the rotated basis is the first 2 lines of Eq. (I3T1) . If we 
linearize the Hamiltonian, then we have following terms 



involving S, 

i i,j 

The first terms describes the rotation of S spin under 
a z direction constant external field for all lattice sites. 
If we only keep this term, we will only get a constant 
dispersion 4J. Since now the ground state is up spin 
for all lattice sites and S 1 can flip up spin to down spin, 
the term S}S} sin 2 9 describes an interaction between the 
flip in site i and neighboring site j . Thus the spin flip can 
propagate and give a plane wave momentum dependence 
like /kS 2 , which agrees with our previous result. The 
third mode corresponds to spin flip: |1, 1) — > \ — 1,-1). 
In order to make this flip propagate, one need a term 
like K^K} 1 in the Hamiltonian. But when J4 = there 
is no such terms, thus the third modes is constant when 
J 4 = 0. When J 4 7^ 0, we have a term K\ x K^ x s\s\ in 
Eq pip , and the third mode will acquire momentum 
dependence. 



2. General Case 

Now we turn to the general case when J4 7^ and 
t' 7^ 0. The t' and J4 terms involve operators like K\\ 
and .K33 which couple the a and r spin and split the two 
degenerate modes we obtained in the decoupled case. We 
denote these two terms by H' . Transforming to the boson 
representation we have 
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H' = [(-2^6 - 2c t c) i sis 2 + (-cU - d f c)iSic 2 + (-6 f d - d f b) l c 1 s 2 + (b f c + c*b) z ac 2 

i 

-2 J 4 ^2 [2(-2& t & - 2^c) l c\c 2 2 + 2(c t d + d) c) lC \s 2 c 2 + 2(bU + d'fy^^cl + 2(b t c + c^b) i sxc l s 2 c 2 

i 

+ Y,(b + b f Ub + b^jcfsl + J2(c + c%{c + J) 3 s\cl + J2(d + J)i{d + d)) jS \sl 
3 3 j 

+2 ^(b + tf)i(c + c r ) j sic 1 s 2 ci - 2 ^(b + tf) % {d + d^jSidsl - 2 J^(c + c)) t {d + d r ) j s 2 1 s 2 c 2 . 
3 3 3 

We can make use of identities like — ts% — t' s\s 2 + AJ\c\ + AJic\c\ = 4(Ji + J4C 2 ,) to simplify the above equation. 
Transforming to momentum space and combining with Hq, the total Hamiltonian is (we drop the subscript k to 
simplify the notation) 

H = ]T( J 2 + Jicl)[8tfb - 2s 2 f k (b + &t) 2 ] + ( J t + J 4C 2)[8ct c - 2s 2 J k (c + c*) 2 ] 

k 

+ [8( Jic 2 + J 2 c 2 ) - 2i(«i + s 2 )l dU - 2J 4 s 2 s 2 f k (d + Sf 
+2(t + AJ lSl ) Cl (bU + d)b) + 2{t + U 2 s 2 )c 2 ( 



-AJ 4 



S1C1S2C1 



Mb + 6t)( c + c t) _ SlClS 2 / k (6 + 6+)(d + d f ) - s 2 s 2 c 2 h{c + c t)(d + d f ) 



(37) 



with /k = (cosfca; + cosfc y )/2. We can rewrite the 
Hamiltonian (|37p in a matrix form H — ip^Mip with 
tp = Qr >b,c ,c,d\d) T and the 6 by 6 symmetric M. 
There are only 12 independent elements because the ma- 
trix elements are the same under the interchange 1 <-> 2, 
3 -H- 4 and 5 ■<-> 6 for both indices or by interchanging 
the position of the two indices ij O ji. These matrix 
elements are easy to read off from Eq. (|37|) . 

Since the Hamiltonian is quadratic, one can use stan- 
dard Bogoliubov transformation to diagonalize the above 
Hamiltonian and find the collective mode dispersions. 
We introduce the following quasi-particle creation and 
annihilation operators 

b = una + vxia' + u\ 2 j3 + V12P + "137 + «137 
c = u 21 a + v 2 \a* + u 22 j3 + v 22 ft + u 2 37 + «237 t 
d = u 31 a + v 3 iaJ + u 32 (3 + v 32 ^ + u 33 ~f + u 3 37 t - 

This introduces 18 real parameters u^^Vij for i,j — 
1, 2, 3. In order to make the quasi-particles to be bosons 
we require that the only non zero commutators are 
[a, a'] — [j3, /3t] = [7, 7^] — 1. This leads to the following 
constraint 

3 

^2(u 2 im -v 2 m ) = 1 

m— 1 
3 

m— 1 
3 

^ ^ijkij^im^jm ^im^jm) — 0- (^^) 
m— 1 



for i,j,k — 1,2,3. There are 9 constraints, thus we are 
left with 9 free parameters. 

Thus there are only 9 independent elements in Hamil- 
tonian so that the off-diagonal 9 can be set to zero by ad- 
justing the 9 free parameters among u and v. Then we are 
left with a diagonal Hamiltonian and 3 dispersion modes. 
But it is very difficult to explicitly solve the 18 real pa- 
rameters u and v. It is easier to determine the the disper- 
sions directly. But this can not be done by directly diago- 
nalize M, because u and v do not form a unitary matrix. 
Introducing <j) = i a \ a , P\ Pi l\ li ) T > tne Bogoliubov 
transformation can be rewritten as tjj — U <fi. One can 
verify that U'U ^ 1. This is because tp'ip ^ <p^4> and ip 
cannot be treat as the norm of the vector tp. This can be 
fixed by introducing matrix a = diag{ 1,— 1,1,— 1,1,— 1}, 
then we have ip'aip = [b, b^] + [c, c 1 } + [d, d''] — 3 and sim- 
ilarly (p^acp = 3. Thus we have U^aU = a. This can 
also be written as <rW ■ all = I and all ■ aU^ = I. One 
can verify that the last equation is consistent with the 
constraints in Eq. (|38jl. 

Suppose we have determined the parameters u and v 
that diagonalize the Hamiltonian, then we have WMU = 
diagjwi, u>i, u> 2 , u> 2 , lu 3i uj 3 } and wi i2 ,3 are the three exci- 
tation modes. We can rewrite this as 

crW ■ Ma ■ crU = diagjwi, — u>\, u> 2 , — uj 2 , uj 3 , — uj 3 } (39) 

It is easy to see that 

det (aU 1 -Ma-aU-ujI) 

= det(cr[/ t )det(Mcr - ul)det(all) = det(Mcr - luI) 

thus matrix aW -Ma-aU and Ma have same eigenvalues. 
We can determine the dispersion of the three modes by 
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solving the equation 

dct(Mcr -ul) = (40) 

which is a cubic equation in u 2 . We determine the roots 
numerically. We again find 3 collective modes but now 
all of them are weakly dispersive. 

Actually we are more interested in the special case 
J\ = Ji- In this case, one can further simplify the equa- 



tion by introducing symmetrized and anti-symmetrized 
variables a\ = (b + c)/ \[2 and a 2 = (b — c) /\/2. Then the 
6 by 6 matrix will decompose into one 4 by 4 matrix for 
symmetrized operator oi and d and one 2 by 2 matrix 
for anti-symmetrized operator a 2 . Now the eigenvalues 
of symmetrized operator a\ and d are 



Wl,3 = 



2{A s A d + B s B d + 2C s C d ) ± 2y/(A s A d + B s B d + 2C s C d ) 2 - 4(A S B S - C 2 )(A d B d - C 



1/2 



(41) 



withA Si<J = Ai±A 2 and the same for B Sid and C Sid . Here we have Ai = 2J(2-s 2 / k )-2J 4 c 2 [-2 + (l + 2/ k ).s 2 ]+t'c 2 /2, 
A 2 = —2(J + 2J 4 c 2 ).s 2 / k , Si = 8Jc 2 - 2ts - 2J 4 s 4 / k , £ 2 = -2J 4 s 4 / k , Ci = V%c(i + 4Js + 2 J 4 s 3 / k ) and C 2 = 
2v / 2J 4 ,s 3 c/ k ; with s = sin# and c = cos#. The eigenvalues of anti-symmetrized operator a 2 is 

1/2 



W 2 



(8 J + 2J 4 c 2 (2 + s 2 ) - t'c 2 ) (8 J + 2J 4 c 2 (2 + s 2 ) - t'c 2 - 8Js 2 / k ) 



(42) 



In order to get a qualitative picture of these three 
modes, we can numerically compute some typical eigen- 
vectors. For example, take the point at zone boundary 
k x = n/a, k y = 0. The three eigen-modes are created by 
the following quasi-particle operators 

a* =0.67(6+ + <:*)- 0.3^, 

ft = 0.71(6+ - c+), 

7+ = 0.21(6+ + c+)+0.95d+ 

One can see that the first and third modes symmetrically 
mix |1,1) -> |1,-1) + |-1,1) and |1,1) — > | — 1, — 1> . The 



first mode has more weight on the former and third one 
has weight on the latter. On the other hand the second 
mode mixes 1 1, 1) — > |1,— 1) — | — 1,1). anti-symmetrically. 

In the decoupled limit, we have two degenerate Ising 
modes for bosons b and c. In other word, these two modes 
are propagation of 90 degree flip. The third mode is con- 
stant corresponding to boson d or 180 degree flip. In 
general case with non-zero J 4 and t', there is no degen- 
eracy and the third mode is also dispersive. In order to 
understand the qualitative effects of J 4 , we first consider 
that t' = and J 4 is small and then make a perturbation 
expansion as follows 



wi = V^J 2 - ht 2 + 



w 2 



= \/l6J 2 - f k t 2 + 



uj 3 = 8 J + 



16 J 2 - 1 2 
\/l6 J 2 - ht 2 
32 J 2 - ht 2 
v/16 J 2 - f k t 2 

A 



1 - (1 + /k) 

1 3t 2 



t 4 



+ 



32J 2 
t 4 



2 64J 2 1024J 4 



1024 J 4 
J 



Ji 
J 



t 2 t 
2J 2 ~ (4 + /k) 128J 4 



We can also assume that J 4 = and t' is small then make a perturbation expansion as follows 



wi = a/16 J 2 - / k t 2 + 



J 2 



w 2 



= V16J 2 - / k i 2 + 



\/l6J 2 - / k t 2 L 
J 2 

\/l6J 2 - At 2 



-2 + (2 + 5/ k ) 



16J 2 



t 4 



2 + (-2 + 3/ k ) 



16J 2 



+ /k 



256J 4 

t 4 
256J 4 



J 
J 



-s = SJ+— 2 t' 



Qualitatively, both J 4 and t' split the dispersions of u)\ these two modes. More precisely, positive J 4 makes u)\ 
and w 2 , but J 4 also changes the relative dispersion of 
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k xy (n/a) 

Figure 5: The collective mode dispersions of Ashkin- Teller 
model with negative J4 as a function of k x along the k x = k y 
direction. The black curve stands for uji and red curve stands 
for CJ2- The parameters are taken to fit the experimental 
results in Hgl201, T c = 65K. The blue curve gives the calcu- 
lated third mode. 



more dispersive than u>2 and negative J4 makes u>\ less 
dispersive than w 2 - On the other hand, t' mostly shifts 
the whole curve and makes both the modes more disper- 
sive. Positive t' pushes uj\ up and push ui2 down, while 
negative t' does the opposite. 



IV. COMPARISON WITH MEASURED 
DISPERSION 

We Choose parameters to reproduce the experiments^ 8 - 
in Hgl2Ql with T c = 61K. The results are shown in Fig. 
([5]). In the experiments, one of the collective modes is at 
40 ± 5 meV another one at 50 ± 5 meV at k = 0. The dis- 
persion across the Brillouin zone in both the (11) and the 
(10) directions for the higher energy mode is 5 ± 5 meV, 
while the lower mode is even less dispersive. We also have 
the very important constraint from the thermodynamics 
of the AT mode l 24 ' 28 that to have a transition with no 
divergence in the specific heat — 1<J4/J<0,J>0. In 
Fig. ([5]), the following parameters are used: 

i = 3.9meV, J = 7.25meV, 

J 4 = -0.2J, t' = -2.15*. (43) 

These give that sin# w —0.257. The dispersion width 
of the first mode is about 10% of the energy gap and 
second one is around 5% of the energy gap. The results 
are shown in figure O These dispersions change with 
doping; we use a doping where a lot of data is available. 

As discussed in Appendix D, the highest energy branch 
with eigenvalue CJ3, is primarily an excitation with angu- 
lar momentum 2 with a small admixture of angular mo- 
mentum 1. The former are undetectable in neutron scat- 
tering experiments. This branch also has a higher energy 
than has currently been addressed by inelastic neutrons; 
it may also be an over-damped excitations since it lies 
well above the pseudo-gap energy in the single-particle 
spectra. 



CONCLUSIONS 



C. Damping 

The principal source of the damping of the collective 
modes is the decay into particle-hole pairs, of similar en- 
ergy and same symmetry as the collective modes, pro- 
vided by Hf-coII- There exists a pseudo-gap in the 
single-particle fermion spectra in the same region of the 
phase diagram as the observed order. Therefore for en- 
ergies of up to about 50 meV in the well under-doped 
region, damping of the collective modes is expected not 
to obscure their visibility. Still, it is surprising that the 
two observed branches of modest 8 - in _ffgl201, T c = 61K 
have a damping smaller than the experimental resolution 
of about 5 meV. This may be due to the large magnitude 
of the order parameter so that most of the orbital cur- 
rent spectral weight of the appropriate symmetry which 
in the non-interacting model is in the incoherent fermion 
excitations is transferred to the collective excitations. It 
is possible that the third branch, which is at about 100 
meV has a significantly higher damping. 



We find that with reasonable parameters, one can fit 
the measured dispersion of the rather unique collective 
modes in the loop ordered phase. The experimental dis- 
covery of these collective modes, whose appearance as 
a function of temperature coincides with the transition 
temperature of the loop ordered phase and whose inten- 
sity is consistent in its variation as the square of the 
measured order parameter as a function of temperature, 
adds further confidence to the existence of loop order. 
Indeed, there is no calculation we know of in which other 
proposed form of order give multiple weakly dispersive 
collective modes. 

The energy of all the modes calculated should — » for 
momentum q — > 0, as T — > T*(x), the loop order tem- 
perature. We suggest experiments to verify this. It also 
follows that the spectral weight of these modes comes at 
the expense of the unusual local quantum fluctuations 26 
in the marginal fermi-liquid region of the phase diagram 
of the cuprates. Accordingly, we expect a strong diminu- 
tion of the low energy part of such fluctuations in the 
pseudo-gap region of the phase diagram. 

When quantum terms are included, the ground state 
is a product over all unit-cells of the linear combinations 
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of the four classical configurations. Beyond mean-field 
calculations, the ground state is more complicated and 
includes effects due to the zero point fluctuations of the 
collective modes. Since the modes are Ising-like, this is 
a small change. We have estimated that the quantum 
correction to the expectations (S 3 ) (T 3 ) and (K 33 ) are 
smaller than 1%. 

The mean-field ground state discovered here has im- 
portant bearing on the effective angle with respect to 
the c-axis, deduced by elastic polarized neutron scatter- 
ing experiments 16 . This is discussed in an accompanying 
paper. It also may have consequences for the gap in one- 
particle fermion spectra in the pseudo-gap regime. We 
hope to discuss this important matter in the near future. 
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same is true for 77 = 171 ■ ■ ■ je — °' z ® I a- The 8 dimen- 
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spinors according to 77 = ±1. The rotation operator can 
be written as 



U = exp 



ft* 



T 1 - K 23 + K 32 — K 11 ) (A3) 



The corresponding operator acting on the spinor space is 



U' = exp 



S 45 + ^25 ~ S 14 ~ S 36 ) (A4) 



So this is a 7r/2 rotation of a 2d plane inside a 6D space. 
Then U is a 2ir rotation in the 6D space. When it act 
on the spinors, it gives an extra minus sign, just as in the 
4D Dirac spinor case. The same thing happened to U 4 
since it act on the chiral spinors. 
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Appendix A: AT model in SO(6) spinor 
representation 

In this appendix, we will rewrite the Ashkin- 
Teller model generators in terms of 50(6) spinor 
representation 27 . In this formalism, one can explic- 
itly show that the spinor representation that requires 
U 4 = — 1 instead of U 4 = 1. From group theory we 
know that locally ST/ (4) is equivalent to SO(6). More 
precisely, SU(4) is a double cover of SO(6). Thus the 
fundamental representation of SU(£) we have used in the 
paper corresponds to the spinor representation of SO(6). 
We can express U in terms of this spinor representation. 
The minimal dimension of gamma matrices of SO (6) is 
8. One can use the eigenvalues of 77 to define chiral 
spinor which is dimension 4, same as the fundamental 
representation of SU (4) . We can introduce the following 
6 gamma matrices 

7a = {a x ® S 2 , cr x ®S 3 , a x <S)K 14 , 

a x ® K 12 , a x ® K 13 , a y ® 7 4 ) (Al) 

and one can verify that they satisfy {j a , 7b} = 2S a b- Then 
the 15 generators of SO(6) are given by S a b — 47 [7a, 7b] 
and they correspond to the 15 generators of 577(4). They 
can be explicitly written as 

Sab = i(/<8> S\ I®T 4 - 2 ' 3 , I®K l \ 

a z ® S 2 ' 3 , o z <g> K 1 ^ (A2) 

with i = 2,3 and j = 1,2,3. Note that all the above 
generators can be decomposed as two 4 by 4 blocks. The 



Appendix B: Need for SU(4) formalism 

In order to treat the 4 states on the equal footing, we 
enlarge the the SU(2) x SU(2) group to SU(4) group. 
Actually, in the classical AT model, only S 3 , T 3 and 
K 33 appear, and all of them are diagonal matrices and 
commute with each other and can be treat as numbers. 
Thus at the classical level there is no difference between 
SU{2) x SU(2) and SU(4) formalism. But for the quan- 
tum model, this is no longer true. In the SU(2) x SU(2) 
formalism, all the quantum states are direct product 
S 2 x S 2 and they can be parametrized by 4 real num- 
bers. There are only two sets of raising and lowering 
operators which can be used as quantum flipping terms. 
In SU{4) formalism, the space of all possible quantum 
states is CP3 which is parametrized by 6 real numbers. 
There 6 sets of raising and lowering operators and they 
exhaust all possible rotations in CP3. 

Therefore SU(2) x SU(2) can not make the most gen- 
eral unitary evolution in the 4 state system. Thus It can 
not treat all the rotations among the 4 states on the same 
footing. In the linear order, the only quantum flipping 
term is a x or t x or equivalent terms, which can switching 
between a z — ±1 or t 2 = ±1. Other types of rotation 
such as cyclic permutation U can not be expressed as 
linear combinations of SU(2) x SU(2) generators. They 
have to be expressed in terms of high order power terms. 

Appendix C: The collective modes in the 
SU(2) x SU{2) Formalism 

Nevertheless, It is interesting to compute the collec- 
tive modes of quantum AT model without enlarging the 
basis to SU(4) but restrict to SU(2) x SU{2) group and 
compare the results. In this way, a 1 and r l are just 2 
by 2 matrices. They correspond to S l and T l matrices 
in the SU (4) basis. There is no operators correspond to 
K l 3 in SU(2) x SU{2) basis and all the quantum states 
are direct product of two sets of Ising spin states. Now 
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Mr/a] 



Figure 6: The two modes in the SU(2) x 517(2) formalism. 

we only have to introduce two HP bosons thus there are 
only two modes. 

Now we introduce the HP boson as follows. 

a x = Uib + tfU u a v = -i(Uib-tfUi), 
a z = 1 - 2tfb 

t x = U 2 c + c f U 2 , r y = -i(U 2 c - JU 2 ), 
t z = 1 - 2c f c 



with U=(l- tfb) 1 / 2 and with J7 2 = (1 - etc) 1 / 2 . 
We consider the simplest case. 

ff<JAT=5>(o?+7f)] 

- E t ^ °i + jT ' T / + J 4<r?a-|r?r?] (CI) 
Now the ground state is 



\G) 



cos-|l) +sin-| - 1) 



cos-|l) +sin-| - i; 



(C2) 



with 9 being determined by t + 4Js + 4 J 4 sc 2 = with 
s = sin 9 and c = cos 9. 

As before, we should first make a basis rotation to 
make the ground state satisfying (a z ) = (t z ) = 1. Then 
we can plug in the HP representation and expand to 
quadratic order. In the momentum space, we find 



H = Yl ( - 2tsb l bk - 1J[-Mb ] k b k )c 2 + (b k + b{) 2 f k s 2 } 

k 

-2Uc\c k - 2J[-4( C t Cfc )c 2 + (c k + 4) 2 /fcS 2 ] 
- 2J 4 [-4c 4 (^ fe + 4c fe ) + 2s 2 c 2 (b k + b\){c k + 4) + (b k + b\ + c k + cl) 2 f k s 2 c 2 ' 

Let i\) = 6, cT, c) T , then the Hamiltonian can be written as H = ip'Mip with 



M = 



Here J' = J + J 4 c 2 



/ 4J'(2- S 2 //c) 
4J's 2 / fc 



-4J's 2 /, 



-4J 4 s 2 c 2 (l + f k ) -4J 4 s 2 c 2 (l + f k ) 



4J'(2-s 2 / fe ) -4J 4S 2 c 2 (l + / fc ) -4J 4S 2 c 2 (l + / fe 



-4J 4S 2 c 2 (l + / fc ) -4J 4S 2 c 2 (l + / fc ) 4J'(2-.s 2 / fe ) -4J's 2 / 
V-4J 4S 2 c 2 (l + / fc ) -4J 4S 2 c 2 (l + / fc ) -U's 2 f k 4J'(2-s 2 / fc ) 



Then we can diagonalize the above matrix and find the dispersion of the two collective modes 



0J k = 8 X /(J + J iC 2 ) (J + J 4 c 2 )(l - s 2 f k ) ± J 4 s 2 c 2 (l + f k ) 

I 



(C3) 



(C4) 



(C5) 



If we let J 4 = 0, then we get oj k = 8J(1 - s 2 / fe ) 1/2 
which is the same as the Ising model solution. Since 
fk = — 1 for k x = k y = it, the two modes are degenerate 
at the zone corner. 

We can still use the parameters from section Mil t 
;').!). J = 7.25, J 4 = -0.2J and t' = 0. Then we find 
the ground state parameter sin6> = —0.34. The two col- 
lective mode dispersions are plotted in figure El In this 
formalism, the two modes are always degenerate at the 
zone corner and the experimental results are not well 



represented. 



Appendix D: Inelastic neutron scattering 

We present here the theory for the inelastic neutron 
scattering from the collective mode of the loop ordered 
phase. For calculating elastic neutron scattering^, We 
approximate the loop current states as 4 spin 1 local mo- 
ments. Based on the same formalism, we can also com- 
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pute the inelastic neutron scattering amplitude. 

First, we need to express the local moment operator in 
terms of the generators of Ashkin- Teller model. We can 
make use of the local moments expression of the 4 loop 



current states^. Then the local spin operator at Ri can 
be obtained by L R — (ip, R4 \L l \ip' , Ri), with ip,ip' take 
the 4 possible loop current states. 



L x Rl = £r 2 = liS 1 + T 1 ), L Rs = L Ri = -liS 1 + T 1 ), 

Lk = i%* = \is 2 + t% L y R2 = l ra = i(5 2 - T% 
l z Ri = -£r 3 = ^ 3 + t3 )> L k 2 = -Liu = \(s 3 - r 3 ) 



Then the inelastic magnetic differential cross section can be expressed in terms of the correlation function of these 
local moment operators 

rPrr V f°° 

jndV = lEE^-^)^'' 11 *^ / (L^ML^e-^dt (Dl) 

Hv a ,b 00 

As in the calculation of the collective mode, it is easy to compute the corelation function in terms of Holstein-Primakoff 
bosons. To the linear order, we find 



L x Rl = £r 2 = ^(6 + & f + c + ct), L Rs = L R4 = -^f(b + 6t + c + c t), 
L v Rl = L R3 = - b+ ct - c), L y R2 = L R4 = - i(6t _ b + c - c t), 

= -On = ~0> + 6* + c + ct), L R2 = -L R4 = + &t _ c _ c t). 



Note that only HP boson operator w and ct appear in the above equation and there is no <P . This is consistent with the 
fact that the local loop current is approximated by spin 1 object. Since S flips S z by 2, it does not couple to neutron 
scattering. We can introduce symmetric and anti-symmetric combinations a% = (6 + c)/y/2 and a2 = (b — c)/v2 to 
simply above equations. The highest energy mode is a mixture of a\ and d\ therefore the highest energy mode should 
still be observable in principle by neutron scattering with reduced amplitude through exciting a\ bosons. 

For notational simplicity, we introduce the following quantity. 

/oo 
{L^ a {Q)L^ b {t))e- lbJt dt (D2) 

a,o 00 

Then the cross-section can be written as 

- [(1 - ql)c* x + (1 - foe™ + (1 - e z )c zz 



dVldE' 1 

-q x q y (C XV + C VX ) ~ %UC yz + C zy ) - q x q z (C zx + C xz ) (D3) 
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We find the diagonal elements of C MI/ as follows, 



C xx = 2 cos 2 9{sm q ■ R4 + sin q ■ R 2 ) 2 



(oi(O)oi(t)) + (a{(0) ai {t)) + (oi(0)o{(t)) + (a\(0)a\{t)) 



C vv = 2 cos 2 q ■ Ri 



+2 cos q ■ R 2 



(ai(O)ai(t)) - (a\{0) ai {t)) - (ai(O)oJ(t)) + (4(0)4(0) 



(a 2 (0)a 2 (*)) - (4(0)« 2 (t)> - (a 2 (0)4W> + (4(0)4(*)> 



'dt 



C zz = 2 sin 2 9\ sin 2 q- Ri 



(ai(O)ai(t)) + (4(0)aiW) + («i(0)4(i)) + (a{(0)a{(t)) 



- sin q • R2 



(a 2 (0)a 2 (t)} + (4(0)a 2 (t)) + (a 2 (0)4(*)) + (4(0)4(*)> 



(D4) 



Here we have made use of the fact that (aia 2 ) = (a|a 2 ) = (ai4) = (44) = 0- 
For the off-diagonal elements, we find that 



Qxy 



C yx oc 



(ai(0) ai (t)) + (ot(O)oi(t)) - (ai(0)4W> - (4(0)4(0) 
(oi(0) ai (t)) - (4(0)aiW> + (ai(0)4W> - (4(0)4(0) 



e^dt 



The cross term of above two equations cancels out. In Eq (|D8|) . we show that (aidi) — (a\a\) and (a 2 a 2 ) 
Because of this identity, the square terms also cancels out, thus we have C xy + C yx = 0. 

Similarly we find that 



- /„t„t\ 



C zy,yz K 



( ai (o) ai (t)) ± (4(o) 0l (t)) t (ai(o)4w> - (4(0)4 w> 

(a 2 (0)a 2 (t)) ± (4(0)a 2 (0) T (« 2 (0)4W> - (4(0)4W> 



e~ lU)t dt 



e~ iu "dt 



Due to the same reason, we have C zy + C yz — 0. Then the only non-zero off-diagonal term is 



C xz = C zx = 2 sin 6» cos 6» (sin q- R x +sinq- R 2 )sin(q- Ri) 

(oi(o)oi(*)) + (4(o)aiW> + (ai(o)4w) + (4(°)4 (*)> 



'dt 



(D5) 



The HP boson can be expressed in terms of Bogoliubov bosons as follows 

ai = uia + V\cv + U37 + «37^ 
a 2 = u 2 (3 + v 2 ft 



(D6) 
(D7) 



with momentum dependent coefficient Ui and vi for i = 1,2,3. Here cv , ft* and 7^ are the creation operators of the 
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three collective modes. At T = 0, we find the following correlation functions 



f 

L 
L 



-iujt 



(ii(ai(O)aiO)) 



dt(a\(0)a\(t)) = uivi6(w - Wi(q)) + u 3 v 3 S(uj - w 3 (q)) 



e~ 4 "*^(ai (0)a\(t)} = u\8{uj - Wi(q)) + u§5(a; - w 3 (q)) 
e" 4 "*^(a t 1 (0)ai(t)) = u 2 £(w - Wi(q)) + uf<J(w - w 3 (q)) 

) 

/•CO 

e- 4 " t ^(a 2 (0)a 2 (t)) = / e- iut dt(4(0)al{t)) = u 2 v 2 S(lu - w 2 (q)) 



./ — ; 



e- 4 "**(a 2 (0)4 (*)) = u 2 2 S(lu - LU 2 (q)) 
e- wt dt{al{Q)a 2 {t)) = v 2 <5(w - w 2 (q)) 



(D8) 



Since we ignored the damping effect, there are only delta function like peaks. Note all 3 collective mode dispersions 
appear in the above correlation functions. Collecting all results, we find 



dVLdE' 



k 



4(1 - q 2 x )F xx + (1 - q 2 z )F zzl - q x q z F x 



(ui + vi) 2 5(lj - wi(q)) + (u 3 + v 3 ) 2 S(lu - W3(q)) 



+ (1 - ql)F yy i (ui - vi) 2 S(oj - wi(q)) + (u 3 - v 3 ) 2 S(oj - w 3 (q)) 
+ [(1 - €)Fyy2(u 2 - v 2 f + (1 - q 2 z )F zz2 (u 2 - v 2 ) 2 }5(uj - w 2 (q)) 



(D9) 



with F xx = 2 cos 2 ^(sinq • Ri + sinq • R 2 ) 2 , F yy i = 2 cos 2 (q • Ri), F yy2 — 2cos 2 (q • R 2 ), F zz \ = 2 sin 2 6'sin 2 (q • Ri), 
F zz2 = 2 sin 2 6*sin 2 (q • R 2 ) and F xz — 2 sin 9 cos 9 (sinq • Ri + sinq • R 2 ) sin(q • Ri). 
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